Ballistic magneto-thermal transport in a Heisenberg spin chain at low temperatures 
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We study ballistic thermal transport in Heisenberg spin chain with nearest-neighbor ferromagnetic 
interactions at low temperatures. Explicit expressions for transmission coefficients are derived for 
thermal transport in a periodic spin chain of arbitrary junction length by a spin-wave model. Our 
analytical results agree very well with the ones from nonequilibrium Green's function (NEGF) 
method. Our study shows that the transmission coefficient oscillates with the frequency of thermal 
wave. Moreover, the thermal transmission shows strong dependence on the intra-chain coupling, 
length of the spin chain, and the external magnetic field. The results demonstrate the possibility of 
manipulating spin wave propagation and magneto-thermal conductance in the spin chain junction 
by adjusting its intra-chain coupling and/or the external magnetic field. 
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I. INTRODUCTION 

Thermal transport properties of low-dimensional sys- 
tems have gained much attention recently. Intensive 
studies in low dimensional systems have made some im- 
portant progress not only in understanding the underly- 
ing physical mechanism but also in controlling of heat. In 
particular, several conceptual thermal devices have been 
proposed such as thermal rectifiers/diodes^, thermal 
transistors^, thermal logical gates^, thermal memory^, 
and some molecular level thermal machines^ 7 -. Much 
work has also been done to the quantum transport in 
nanostructures^. 

Low-dimensional systems, especially one-dimensional 
(ID) materials, offer the possibility to study quantum 
effects that are masked in three-dimensional systems. 
In recent years, many interesting experiments on ther- 
mal transport in ID spin chains 9 - ! 10 ' 11 ! 12 are performed, 
where the lD-spin-chain compound materials give us nice 
physical realizations of ID toy 1 model systems. From 
these experiments, it is possible to control heat trans- 
port in spin systems by a magnetic field. There are also 
theoretical studies about thermal transport in ID spin 
chains, some of which show anomalous transport due to 
the integrabilit y 13 ' 14 i 15 ' 16 , such as anisotropic Heisenberg 
5=1/2 model, the t-V model and the XY spin chain. 
The properties of energy transport through the chains 
differ for different anisotropy of the intra-chain spin in- 
teractions. In all spin systems, the mean free path of 
itinerant spin excitations increases as temperature de- 
creases. Therefore in the low temperature limit, the ther- 
mal transport in spin chain can be regarded as ballistic. 

The phononic transmission coefficients in quasi-one- 
dimensional atomic models can be calculated by transfer 
matrix metho d 17 ! 18 ' 19 ! 20 . However, if there are evanes- 
cent modes with large |A| (A is an eigenvalue of the trans- 
fer matrix), the evaluation of the transfer matrix can be 
numerically rather unstable, particularly when the sys- 
tem size becomes large. Alternatively, nonequilibrium 



Green's function (NEGF) is an efficient method to calcu- 
late the transmission coefficient. Unfortunately, both of 
these two methods can not give an analytical expression 
easily. 

In this paper we give an explicit analytical expression 
of transmission coefficient through a spin-wave model 
from a wave scattering picture. Here we study magneto- 
thermal transport in an isotropic Heisenberg spin junc- 
tion coupled to two semi-infinite spin chains in equilib- 
rium at different temperatures. By Holstein-Primakoff 
transformation we map the spin operators to spinless bo- 
son operators, and consider only the harmonic terms of 
Hamiltonian in the low temperature limit, which is dis- 
cussed in Sec. |TTJ The analytical solution from the spin- 
wave model is shown in Sec. IIII1 we get the explicit for- 
mula for transmission coefficient. In Sec. IIV1 we intro- 
duce the nonequilibrium Green's function method, and 
use it to study thermal transport. The results and dis- 
cussion are given in Sec. El A short summary is presented 
in Sec. ED 

II. MODEL 

The Heisenberg spin chain consists of three parts: two 
semi-infinite leads and an arbitrary junction region (see 
Fig. [1]). The two leads are in equilibrium at different 
temperature Tl and Tr. We apply different external 
magnetic field to the three parts along z direction. So 
the Hamiltonian of this system is given by 

i i 

here Ji is the interaction between spin site i and i + 1, 
hi is the magnetic field applied to spin site i. Using 
Holstein-Primakoff transformation 2 ^ 

S + = V2S—a + a a; S~ = a + \/2S— a+a; S z — S — a + a, 

(2) 



2 



it is easy to map spin operators to spinless boson opera- 
tors a + , a. In the low temperature limit, (a + a) <C 2S 1 , we 
obtain the following Hamiltonian by neglecting the terms 
containing products of four or more operators, 



(3) 



where £ = ~E ( J i s2 + h i s )> and K u = (A-i + Ji)S + 

i 

hi, Ki t i+i — Ki + \_i = —JiS, here S, which can be any 
integer and half-odd-integer, is the maximum value of 
spin. We choose 5=1, without loss of generality. 
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FIG. 1: (color online) The system is an infinite Heisenberg 
spin chain, which consists of two semi-infinite leads with an 
arbitrary junction region. The two leads are in equilibrium 
at different temperature Tl and Tr. We can apply different 
magnetic fields to the three parts. 



III. ANALYTICAL SOLUTION FROM THE 
SPIN- WAVE MODEL 

Because only harmonic terms are contained in the 
Hamiltonian, we can assume a spin wave solution trans- 
mitting from the left lead to the right lead through the 
junction region. We consider the two leads as uniform 
spin chains of intra-chain spin coupling Jl in a magnetic 
field h^. The junction has an alternating coupling, Ji 
and J 2 in a field h. The unit cell of the junction part 
contains two spin sites. We assume an incident wave as 



When it arrives at the center part, partially 



will be reflected and partially will be transmitted. The 



reflected wave is rX l e 
can be written as 



and the transmission wave 



iXi 



— 1 „— iut 



vX 



j—l—iut 



(4) 
(5) 



The coefficients r, u, v are obtained from the continuity 
condition at the interface. The transmission wave will be 
reflected and transmitted by the right end of the junction, 
the amplitude of reflection wave are ur u \ and vr v \ for the 
odd site and even one, respectively. The transmission 
wave from the right end is 



iX^hXlt 



(6) 



where L = N— 1, N is the number of sites in the junction 
part. The reflection wave will be reflected (r u2 x 100% 



percent) at the left end again, and then reflected again. 
Finally the total wave function transmitted from the 
junction is a superposition of multiple reflections and 
transmissions, 

ipj = uAf (1 + r ul r u2 X 2 2 L + (r ul r u2 X\ L f + ■ ■ -)Mje- M *. 

(7) 

From time dependent Schrodinger equation 



at 



(8) 



here we set h — 1 for simplicity, we can get the dispersion 
relations for the leads Ai = e lqi and for the junction 

A 2 = e^ 2 as 

1 



uj-{2J L + h L ) = -J L ^—+Xij; (9) 
(w - ( Ji + J 2 + h)f = Jl + 4 + Ji J 2 + xfj . (10) 

Which root of the equations should we use? By adding a 
small imaginary part to lu, that is, replacing it by u + irj, 
then none of the eigenvalues A will have modulus exactly 
1. Considering r\ as a small perturbation, we find for the 
traveling waves 22 



|A| = 1 



a 
v 



n 



(ii) 



That is, the forward moving waves with group velocity 
v > have |A| < 1. In the formulas below, we take 
the root with |A| < 1. The energy band for our model 
(Ji < J 2 ) is 

([Mi + 2Ji)] U [/ l + 2J 2 ,/ l + 2(Ji + J 2 )])n[/i L , h L + AJ L ]. 

(12) 

Finally, we obtain the transmission coefficient (for N 
odd) as 

utiX^ 



r u ir U 2X\ L 



Here, 



- ( Ji + J 2 + h) _ Ji/A 2 + J 2 A 2 
JiA 2 + J 2 /A 2 uj - ( Ji + J 2 

J L (i-x\) 



h) 



TU) '—^r ■ (13) 

(14) 
(15) 

rul = {J L + J-2 + h) + J 2 X 2 /a + J L X 1 ' (16) 

= ^ - ( Jl + Ji + h) + Ji/(aA 2 ) + J L Ai 
'"" w-(J L + Jx + h) + JiX 2 a + J L Xi ' 1 ' 

(18) 

h = Ai(l + r„i), £ 2 = Ai(l + r u2 ). (19) 

If the number of the sites is even, that is, the length of 
the chain is odd, the transmission can be written as 



.—— , v — ua 
Jl - J2 - Ai Jl + aJ 2 /X 2 

lu - (J L + J 2 + h) + aJ 2 /X 2 + J L X 1 



f m = 



uX%t\ 



' ul 



a 2 r u2 



r' ul r u2 Xl L ' 
i'l = at 2 . 



(20) 
(21) 
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If J\ = J2 , all the formulae are reduced to those of the 
uniform spin chain. Although we only discuss the period- 
two spin chain, it is easy to derive similar formulae of 
transmission coefficient for any other arbitrary periodic 
junction by this method. 

For ballistic transport, the thermal current can be writ- 
ten as a Landauer-type expression: 



{!) = 



1 

2^ 



w [f L (u>) - f R (u>)]f(cj)du. (22) 



The conductance is 



— / duj to 1 (iv ) - 



2tt 



dT 



(23) 
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FIG. 2: (color online) The transmission coefficient of the uni- 
form spin chain with coupling J = 1. We apply a magnetic 
field to the junction. The transmission T(ui) shifts along the 
frequency axis with magnetic field and oscillates with fre- 
quency in the range of [1,4]. The number of sites in the junc- 
tion part is TV = 9, with magnetic field h = 1. The scattered 
crosses and solid line correspond to the results from NEGF 
and the spin-wave model, respectively. The inset shows the 
case with a weak magnetic field, h = 0.05. The numbers of 
peaks is TV — 1. 



IV. NONEQUILIBRIUM GREEN'S FUNCTION 
METHOD 

From the discussion in Sec.|IIl we can write the Hamil- 
tonian, by neglecting the ground state energy E$, as fol- 
lows 



STr^+^LC c 

m rn 



' l m~^mt a l) 



h.c. 



(24) 

where H a = £ a" + IfX. a = L,C,R, here 'L,C,R 

lm 

denote left lead, center part and right lead, respectively. 



The Hamiltonian matrix of the full linear system is 

/ K L V lc \ 
H = V CL K c Vcb. ■ (25) 
\ V RC K r ) 

We use nonequilibrium Green's function method^ to 
study the thermal transport in the spin chain. First we 
define the retarded Green's function as 



G r (t,t r ) = -^0(<-^')([a(^) J a + (^')]>• 



(26) 



In nonequilibrium steady states, the Green's function is 
time-translationally invariant and so it depends only on 
the difference in time. The Fourier transform of G r (t — 
t') = G r {t,t') is defined as 



G r [u>] 



+ OO 



G r {t)e iujt dt. 



We also need the advanced Green's function 
G a {t,t') = i6(t' - t){[a(t),a+(t')}), 
the 'greater than' Green's function 

G>(t,t') = -i{a(t)a+(t')), 
and the 'less than' Green's function 

G < (t,t') = -i(a + (t')a(t)}. 



(27) 



(28) 



(29) 



(30) 



Without interaction, the free Green's functions for three 
parts in equilibrium can be written as: 

((ui + ir))- K a )g r a (u) = I, a = L,C,R, . , 

And there is an additional equation relating g r and g < : 

5<H=/HL 9 >)- 5 >)], (32) 

where f(u>) = (a + a) = [e"/ T — 1] _1 is the Bose-Einstein 
distribution function at temperature T; we have set the 
Boltzmann constant Ub = 1- 

The energy current flow from the left lead to the central 
region is 

II = -(Hl), (33) 
which can be expressed by Green's function as 
I L = \im_ Tv{V CL K L G> c (t,t') - G> L (t,t')K L V LC }. 

In frequency domain, the current expression can be writ- 
ten as 

/CO J 
Tr{G> L (u)K L V LC ~ V C lK l G> c (cj)}^. 

(35) 
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FIG. 3: (color online) The transmission coefficient of the junc- 
tion part with periodic coupling Ji = 1 and J2 = 0.5. There is 
no magnetic field applied to the spin chain. The transmission 
oscillates in the band of [0,1] and [2,3], which is consistent 
with Eq. (|12p . Here, the number of sites in the junction part 
is N = 33. The crosses and solid line curve correspond to the 
results from NEGF and spin-wave model, respectively. The 
inset shows the fine curves in the range of ui £ [0.5, 0.6] and 
[2.5,2.6]. The results from two methods are consistent with 
each other. 



FIG. 4: (color online) The transmission coefficient of the 
periodic intra-chain coupling junction in an external mag- 
netic field h — 1. Here, the coupling of leads Jl = 1, and 
Ji = 1, J2 = 0.5 for the junction. The upper and lower panel 
are for N = 33 and N = 34, respectively. The results of the 
spin-wave method (solid line curve) (for N = 33 and TV = 34, 
we use Eq. (|13[) and Eq. (|20p . respectively) is identical with the 
numerical results from NEGF (the scattered crosses) 



Because of the following relations 



(36) 



Gcc = 9c + gc^Gcc, (37) 

S = Si + £r, S Q = VcugaVaC, (38) 

Gcl = GccVcLgL, Glc — QlVlcGcc, (39) 

the current (I) = — (Ir)) can be reduced to 

Landauer-type expression Eq. I|22[). where the transmis- 
sion coefficient is 



f(u) = Ti{G r cc r R G a cc r L }. 



(40) 



The T Q functions are given by T a = — 
For the ID spin junction coupled to two semi-infinite 
leads, which are uniform spin chains of intra-chain in- 
teraction Jl in a magnetic field , the transmission co- 
efficient can be written as 



f(u) = 4Jt(lm(g r 00 )) 2 \G N - 1 , f 



(41) 



where gg = — Ai/Jj,, Ai is given by Eq. ©, Green's 



function G r cc is abbreviated as G and 



Gjv-i,o — (w— K c ~^) N 1 _ l o 



det (oj - K c - S) ' 



(42) 



For N = 3, we can get 6*2,0 = 



Jl fi „ r3 ) here a = 



abc—aj~—cj? 



lu — Jl — Ji — h — J 



2 a r 



LU — Ji — J 2 — h, C = LU 



J^ — J 2 — h — J\g§- For general N, it is difficult to get 
an explicit formula. 



V. RESULTS AND DISCUSSIONS 

In our calculations, we take, ks = 1, h = 1. The unit 
of coupling J is lmeV, then the unit of magnetic field is 
17. 5T; the unit of temperature is 11.6 K; and the unit of 
conductance is 3.86 x 10 _2 nW/K. 

If the whole system is uniform, i.e. the magnetic fields 
applied to the three parts - two leads and the junction 
- are the same, the transmission coefficient T(u>) is al- 
ways 1 in the whole domain [h, h + 4J]. However, if the 
magnetic fields in three parts are different, e.g., h > in 
the junction and = in the two leads, the transmis- 
sion coefficient oscillates with frequency lu in the domain 
[h, 4 J], where the whole system has the same intra-chain 
coupling J. If the magnetic field is weak, the oscilla- 
tion region is very near to 1, and the number of peaks 
(T(lu) = 1) is equal to N — 1. With the strengthening 
of magnetic field in the junction, the transmission T(u>) 
shifts along the axis of frequency lu and the oscillation 
range extends to the domain [0, 1], some peaks will be 
cut off because the shift of curve. The numerical results 
come from NEGF are exactly the same with analytical 
solution from the spin-wave model, which is shown in 
Fig. [2] All the phenomena are still the same when the 
size of the spin chain is very large. However, for small 
size spin chain, the transmission at the forbidden band is 
not zero because of quantum tunneling effect, which can 
be given from both of the two methods and the results 
are exactly consistent. The shift and oscillation of trans- 
mission coefficient is because of interference of the spin 
waves transmission through the junction. 
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FIG. 5: (color online) The transmission coefficient for junc- 
tions connected in series. Here, every segment has 5 sites and 
all segments in a magnetic field of h = 0.2 are connected by 
a spin with no magnetic field, (a), (b), (c) and (d) corre- 
spond to 1,2,8 and 32 segments, respectively. For 32 or more 
segments, there are obvious 5 gaps (forbidden bands ), i.e., 
6 transmitted bands,which is due to the period of magnetic 
field along the chain, here the period is equal to 6 (5 sites 
with magnetic field plus 1 spin with no magnetic field). 



If the intra-chain spin coupling of junction is different 
from that of leads, the transmission coefficient oscillates 
also with the frequency. In Fig. [31 we show the transmis- 
sion coefficient for the junction with periodic coupling: 
J\ = 1, J% = 0.5 and the leads with J = 1, oscillates 
with frequency in the energy band [0, 1] and [2,3] which 
is consistent to Eq. (fT2")) . From this figure, it can be fur- 
ther concluded that our analytical results from the spin- 
wave model is exact. The numerical results from NEGF 
method is consistent with this analytical approach. 

When an external magnetic field ft, = 1 is applied to 
the junction part, the transmission coefficient shifts to 
[1, 2] (J[3, 4], and the oscillation shape changes, which is 
shown in FigfJ] The shapes of oscillation are different for 
odd-site and even-site junction. For the even-site junc- 
tion, the spin chain is symmetric along the chain direc- 
tion, while it is asymmetric for odd-site case. Therefore, 
the reflection in two ends of the junction for asymmet- 
ric chain are different, which causes the difference of the 
transmission compared with symmetric case. 

From the above results, we know that the periodicity 
of the junction can give rise to gaps in the transmis- 
sion. Can we merely apply magnetic field periodically 
to the junction to induce gaps in the transmission, while 
the whole system have the same coupling J = 1, there- 
fore, we can choose the frequency to transmit from the 
junction? It is possible, because that the transmission 
oscillates with frequency when the junction is applied a 
magnetic field. If we connect many junctions in series, 
then the range of oscillation extends and may give gaps 
in transmission. In Fig. at first we let the heat transfer 
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FIG. 6: (color online) The transmission coefficient for a fixed 
junction in different periods of magnetic field. The total num- 
ber of sites is N = 191. The magnetic field is applied to the 
junction as h(p * m) = 0, m = 1,2, 3..., others = 0.5, p is the 
period of magnetic field along the junction, (a), (b), (c) and 
(d) correspond to uniform, p — 2, p = 6 and p — 12 magnetic 
field, respectively. It is shown that there are p — 1 gaps for 
different cases. 



through a 5-site junction in a magnetic field h = 0.2, the 
transmission oscillates a little near 1; if we connect two 
junctions together by a site without magnetic field, that 
is, N = 11, h(i = 1 ~ 5; 7 ~ 11) = 0.2, h(6) = 0, the os- 
cillation will be extended. The gaps are shown evidently 
when 32 or more segments are connected in series by the 
sites without magnetic field. For a fixed size junction, if 
we apply magnetic field periodically, there are p—1 gaps 
( p is the period of magnetic field along the junction) in 
the transmission, which arc shown in Fig. [6l Therefore, 
we can choose the frequencies to transmit through the 
junction by adjusting the periodicity of the intra-chain 
coupling or the magnetic field. 

We can calculate the thermal conductance from trans- 
mission coefficient. In Fig. [7j we show the thermal con- 
ductance versus system size, temperature and magnetic 
field for uniform and periodic spin chain. It is shown that 
the conductance is independent of the system size be- 
cause of ballistic transport while the oscillation of trans- 
mission coefficient changes with the system size. And the 
conductance increases to a constant with the increase of 
temperature, which is consistent with all other ballistic 
cases. From Fig. [Jjb), we see that the conductance de- 
creases to zero as the intensity of magnetic field is in- 
creased. These results indicate that the transmission co- 
efficient will shift along the axis of frequency, which cuts 
off the contribution of low frequency to the thermal cur- 
rent, therefore the thermal current decreases. Because 
the thermal wave of low frequency contributes most to 
the heat flux, the thermal conductance decreases quickly 
with the increase of the intensity of magnetic field. Al- 
though the transmission has a big difference for uniform 
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FIG. 7: (color online) (a) The thermal conductance, a, versus 
system size, N, and temperature T. a increases with tem- 
perature to a constant whereas it keeps invariant with system 
size, where magnetic field is h = 0.5. (b) The thermal conduc- 
tance a versus magnetic field, h and temperature. It shows 
that the conductance decreases to zero with strengthening 
magnetic field at different temperature. The upper blue sur- 
face and lower black surface correspond to uniform (J = 1) 
and periodic (Ji = 1, J 2 = 0.5) spin chains, respectively. 

and periodic chains, the thermal conductance has the 
similar behavior merely with a difference of magnitude. 

VI. CONCLUSION 

In this paper, we have studied the ballistic magneto- 
thermal transport in a Heisenberg spin chain at low tem- 
peratures. We have obtained explicitly an analytical ex- 
pression for the transmission coefficient through a spin- 
wave model from a wave scattering picture. The ana- 
lytical results have been verified by the noncquilibrium 
Green's function method. We have found that the trans- 
mission coefficient oscillates with the frequency because 



of interference of the transmission waves; furthermore, 
the thermal transmission coefficient shows strong depen- 
dence not only on the intra-chain coupling and length 
of the spin chain but also on the external magnetic field. 
There are gaps in the transmission, the number of gaps is 
equal to p— 1, where p is the period of intra-chain coupling 
or the external transverse magnetic field applied to the 
junction, i.e., the number of transmitted bands is equal 
to the value of periods. Therefore, it is easy to choose 
special frequencies to transmit through the spin chain 
junction by adjusting its intra-chain coupling or the ex- 
ternal magnetic field and the heat current in the junction 
can be switched off with the magnetic field strengthen- 
ing. The thermal conductance of Heisenberg spin chain 
at low temperature tends to a constant with the temper- 
ature increasing, decreases to zero with intensity of the 
magnetic field, while it has no dependence on the system 
size. 

Our analytical spin- wave model solution can be applied 
to the ballistic magneto-thermal transport in an arbi- 
trary periodic spin chain. The properties of the magneto- 
thermal transport found in this paper provide the possi- 
bility to manipulate magneto-thermal conductance and 
the propagation of spin waves in the Heisenberg spin 
chain, which may have potential applications in thermal 
control and designing of filter and waveguide for spin 
waves. 
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